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ABSTRACT 

We present spectroscopic analysis of the broad absorption line outflow in 
quasar SDSS J1512+1119. In particular, we focus our attention on a kinematic 
component in which we identify Pv and Siv/Siv* absorption troughs. The 
shape of the unblended phosphorus doublet troughs and the three Siv/Siv* 
troughs allow us to obtain reliable column density measurements for these two 
ions. Photoionization modelling using these column densities and those of He I* 
constrain the abundance of phosphorus to the range of 0.5-4 times the solar 
value. The total column density, ionization parameter and metalicity inferred 
from the P V and S IV column densities leads to large optical depth values for 
the common transition observed in BAL outflows. We show that the true Civ 
optical depth, is ~ 1000 times greater in the core of the absorption profile than 
the value deduced from it's apparent optical depth. 

Subject headings: galaxies: quasars — galaxies: individual (SDSS J1512+1119) - 
line: formation — quasars: absorption lines 
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INTRODUCTION 



Active galactic Nuclei (AGN) outflows have been detected a s blueshifted broad 



absorption line s (BALs) in t 



Dai et al. 



2008 



le UV spectra of ~ 20% of quasars (IHewett fc Foltz 



Kni gge et al. 



Seyfert galaxies ( 



Crenshaw et al 



2003 



2008) an d as narrow absor ption lines (NALs) in ~ 50 % of 



2003 



Dunn et al. 



20081 ). There is growing evidence that 



these ubiquitous sub-relativistic ionized outflows play an important role on sub-parsec as 
well as kilo-parsec scales in controling the growth of the central black hole, the evolution 



of the host galaxy anc 



Elvis 



2006 



Moe et al. 



the c hemical enrichment of the Intergatlactic medium (IGM) (e.g. 



2009 



Ostriker et al. 



20101 ). Moreover, the study of the abundances 



of metals in these outflows observed over a range of redshifts (up to z ~ 5) provide us with 
a unique probe to investigate the history and evolution of the chemical enrichment over 



cosmol ogical sea 



es, which constrains star formation s cenarios and evolution of the host 



galaxy ( 


Hamann 


199 r 


7 . 1998 


: Hamann Sz Ferland 


1999: 


Hamann e 


t al. 


2003: 


Di Matteo et al. 


2004; 


Hamann et al. 


2007 


; ( 


Germain et al. 


2009 




Barai et al. 


I2OII 


)■ 



Studying absorption lines from AGN outflows is the most direct way to determine 
chemical abundances in the AGN environment. This is done by comparing the column 
densities associated with ionized species of the different elements observed across the 
spectrum, combined with photoionization analysis. The major advantage of using 
absorption lines over emission lines in abundance studies resides in the fa ct that they 



19981 ). 



provide diagnostics that largely do not depend on temperature and density (lHamann 
Early abundance studies in BAL outflows implied particularly high abundances of heavy 
elements relative to hydrogen. In several cases, enhancement of carbon, nitrogen, oxygen 
and silicon b y factors of ten s to hundreds of times the solar value s were reported in several 



objects (e.g. 



Turnshek 



1986 



Turnshek et al. 



1996 



Hamann 



19981 ). in contrast to the order 



1 Based on observations collected at the European Southern Observatory, Chile, PID:87.B- 
0229. 
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of magnitude or less, generally derived from the analysis of t 


le quasar emission lines 


Hamann & Ferland 


1993: 


Hamann et al. 


2002; 


Dietrich et al. 


2003; 


Juarez et al. 


2009) 



Perha ps the most puzz 



with P v (IJunkkarinen et al. 



ing_o 



Dservation was the detection of BALs associated 



1995 



Arav et al. 



2001b 



Hamann 



1998 



Hamann et al. 



2003 ). Phosphorus i s ~ 900 times less abundant than carbon in the solar photosphere 



( ILodders et al. 



20091 ). Since Pv and ClV have similar ionization potentials they are 
formed in similar environments. This suggests, based on direct comparison of the measured 
column de nsities, an overabundanc e of phosphorus over carbon of ^ 100 tim es the sola r 



value (e.g. 



Junkkarinen et al 



1995 



1997 



Turnshek et al 



1996 



Hamann 



19981 ) 



Shields 



(119961 ) suggested a scenario consistent with the reported phosphorus overabundances in 



which the enrichment of the BAL material is main 



novae. H owever, our group fjAray 



2001b 



1997 



a 



2002 



2003 



Hamann et al 



Scott et al. 



19971 



1997 



2004 



Telfer et al. 



Arav et al 



Gabel et al 



1998 



y cont ro 



1999bl a: 



led by a popu 



de Kool et 



al 



200 5bf) and others (IBarlow fc Sargent 



Churchill et al. 



1999 



ation of galactic 



2001 



Arav et al. 



Ganguly et al. 



19991 ) 



showed that column densities derived from the apparent optical depth analysi s of BAL 



troug 



is are unreliable due to non-b 



(I19981 ): lHamann fc Ferland! (119991 ): 



ack saturation in the troug 



is. Therefore, 



Hamann et al. 



(12003 ): and 



Leighly et al. 



Hamann 



(I2OO9J) 



suggested that the extreme overabundance of phosphorous relative to carbon is an artifact 
of very high levels of saturation in the C IV troughs, compared to only mild (or non) 
saturation in the Pv troughs. Subsequent measurements of abundances in Seyfert and 
quasar outflows accounted for non-black satu r ation and yielded abundances of only a few 



times solar in the outflows (IGabel et al. 



2006 



Arav et al. 



20071). 



The non-black saturation hypothesis was largely accepted by the community to explain 
the Civ/Pv BAL observations. But this scenario implies that the actual optical depth in 
the C IV trough is roughly 1000 times larger than the apparent one, an assertion that was 
never verified empirically. In this paper, we study the UV outflow of SDSS J1512+1119, 
which exhibits deep absorption troughs from P V as well as S IV. In particular, we report the 
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detection of the excited SlV A1073.51 line, a transition ten times weaker than the excited 
S IV A1072.96. Together with estimates of the number density provided by the analysis 
of absorption troughs from excited states of C ill and Fe III, we pinpoint the S IV column 
density. Photoionization modeling, using the derived column densities as input, shows that 
the phosphorus abundance is close to the solar values. This allows us to confirm that the 
true C IV optical depth is ~ 1000 times greater in the core of the absorption profile than 
the value deduced from apparent optical depth measurements. 

The plan of the paper is as follows: In §[2] we present the VLT/X-Shooter observations 
of SDSS J1512+1119 along with the reduction of the data. In § [3] we identify the spectral 
features and estimate the column density associated with each ionic species. We discuss 
the photoionization solution for the absorber and the implied phosphorus abundance in § HI 
We conclude the paper by summarizing the key points of the analysis in § [5j 



2. Observation and data reduction 



SDSS J1512+1 



Hewett fc Wil 



19 (J2000: RA=15 12 49.29; dec=+ll 19 29.36; z=2.1062 ±0.0020 



2010l : V = 17.7), also identified as Q1510+115 is one of the objects originally 



discovered in the spectroscopic surv ey conducted 



with the UK Schmidt telescope (cf. 



Sargent et al. 



jy Ha zard using objective prism plates 



19881 ). Later spectroscopic observations 



with the double spectrograph at the Palomar Hale telescope revealed the presence of broad 
(Av ~ 1700 km s _1 in Civ) absorption troughs ass ociated with Lyq, C iv and Nv while 



several resolved components were identified in Siiv (ISargent et al. 



1988|). 



We observed the quasar SDSS J1512+1119 with the VLT X-Shooter spectrograph 
on April 26 2011 as part of our program 87.B-0229 (PI: Benn). X-Shooter is the second 
generation, wide band (3000 A to 24000 A), medium resolution (R ~ 6000) spectrograph 
installed at the Cassegrain focus of VLT/UT2. In this instrument the incoming light is 
split into three independent arms, each arm consisting of a prism-cross-dispersed echelle 
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spectrograph optimized for the UV-blue, visible and near-IR wavelengths (UVB, VIS and 
NIR, respectively) which allows for coverage of the full bandwidth in a sin gle exposure 



A detailed description of the instrument and performance can be found in 



Yernetetal 



fl201lh . The total integration time for the UVB, VIS and NIR arms are 8400, 8400, and 
8700 s, respectively. The observations were performed in the slit nodding mode with two 
positions using a slit width of 0.8" in the UVB and 0.9" in the VIS and NIR leading to 
respective resolving powers R = AA/A of 6200, 8800 and 6100. Except for a line from the 
metastable 2 3 S excited state of He I (He I* A3889.80, discussed in a forthcoming paper), no 
additional diagnostic lines are observed within the NIR range so that we limit the current 
study to the UVB and VIS range of the data. 



The observations were reduced in nodding mode using the ESO Reflex 



(Ba llester et al. 



20111 ) and the ESO X-Shooter pipeline version 1.4.5 ( jModigliani et al. 



workflow 



20101 ) in order to obtain the rectified and wavelength calibrated two-dimensional spectra for 
each arm. After manually flagging the remaining cosmic ray hits in the individual frames, 
we extracted one- d imensional sp ectra using an optimal extraction algorithm based on the 



method outlined in iHornd ( 1l986l ). An identical treatment was performed on the observation 
of the spectroscopic standard star LTT7987[ 3 J observed on the same day as the quasar, 
allowing us to flux calibrate the spectra in the range 3200 - 9500 A. We present the reduced 
X-shooter spectra in Figure [H 



3. Analysis of the absorber 

Comparing the unblended line profile of the Siiv doublet, we distinguish five main 
kinematic components associated with the intrinsic outflow in our VLT X-Shooter data. 



2 Reflex is available at |http://www. eso.org/sci/soft ware/pipelines/| 

3 The flux calibrated reference spectrum can be found at: 



http: / /www.eso.org/sci/observing/tools/standards/spectra/ltt7987.html 
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Fig. 1.— Reduced UVB+VIS X-shooter spectra of the quasar SDSS J1512+1119. We 
indicate the position of most of the absorption lines associated with the intrinsic outflow 
(see Section [3]). The VIS spectrum has not been corrected for atmospheric absorption (the 
"+" indicate the location of the major O2 and H 2 atmospheric bands within the spectral 
range). This does not affect our study since our diagnostic lines are located in regions free 
of such contamination. 
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Fig. 2. — Intrinsic absorption troughs associated with Siiv in the VLT X-Shooter spectrum 
of SDSS J1512+1119. At the resolution of the instrument, the outflow is resolved in five 
unblended components spanning a range of velocities between -2100 and -520 km s -1 . The 
continuous line represents the unabsorbed emission model we built (see Section l37Ll) . and the 
dashed lines represents the continuum and emission line models. 
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The centroid of the absorption features are located at radial velocities v ~ -2100, -1850, 
-1500, -1050, and -520 km s" 1 in the rest frame of SDSS J1512+1119 (see Figure |2j). 
Using the Siiv components as a template, we identify absorption troughs associated with 
the kinematic components in a series of ions spanning a range of ionization. While the 
kinematic components are deep and blended together in the ubiquitous Lya, Ly/3, O VI, Nv 
and C IV transitions, some of the components are resolved in the Mg II, C II, Si II, Si ill, Al II, 
Aim and He I* lines (see Figure [3]). We report the detection of deep absorption troughs 
associated with the P v doublet in the component located at -1850 km s _1 (component 
2, see Figure E]), as well as absorption associated with ground state and excited Siv, and 
excited C ill and Fe ill. Of particular interest in this component is the detection of the high 
ionization excited S IV A1073.51 transition, a line ten times weaker than the S IV A1072.96 
arising from the same excited level. The combined detection of these lines allows us to 
accurately determine the total Siv column density for that system (see Section [372]) . In the 
remainder of this analysis we will focus our attention on component 2, deferring the study 
of the other components to a forthcoming paper. 



3.1. Ionic column densities 

In order to compute the column densities, we first need to estimate the spectrum of the 
background source (continuum and emission lines) illu minating the ab s orbin g material. We 



model the unabsorbed and dereddened (E(B-V)=0.051 ISchlegel et al.l ( 119981 )) continuum 
by a single power law of the form F x = Fn o(A/1100) a , where Fi 1Q0 = 77.01 ± 0.05 x 10 17 
erg/cm 2 /s/A is the observed flux at 1100 A (rest frame) and the power law index 
a = -0.403 ± 0.004. Given the overall high S/N of the X-Shooter data (S/N ~ 30 - 70 
over most of the UVB/VIS range), we then fit the unabsorbed line emission using a smooth 
third order spline fit. We constrain the emission model over the broader troughs associated 
with Ly/3, O VI, Lya, Nv, and Civ by using the profiles derived from the unabsorbed red 
wings of each emission line. This fitting procedure is particularly effective given the tight 
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Fig. 3. — Normalized absorption line profiles of the absorption troughs associated with the 
UV outflow in SDSS J1512+1119 for kinematic component 2. For unblended doublets, we 
overplot the expected residual intensity in the strongest component based on the residual 
flux observed in the weakest component assuming the absorbing material completely and 
homogeneously covers the emission source {lotted = Ir, where Ir is the residual intensity in 
the weaker component of the doublet). 
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constraints on the unabsorbed emission due to the high S/N and the presence of mainly 
narrow absorption features across the spectrum. 

Once normalized by the emission profile, we extract the ionic column densities AT [ion] 
for each species by modeling the observed residual intensity Ii (see Figure |3]) inside the 
troughs as a function of the radial velocity We use three different models in order to 
account for inhomogeneities in the absorbing material; apparent optical depth (AOD), 



parti al covering (PC) and a power law model (PL) (see 



Edmonds et al 



2011 



2012 



2004 



or details). We integ rate the line of sight averaged column density (lHamann fc Sabra 



Edmonds et al. 



Borguet et al. 



20111 ) over the radial velocity range v G [—2030, —1570] km s _1 
corresponding to kinematic component 2 and report the values in Table [TJ We only report 
lower limits derived with the AOD model on column densities for singlet lines (Al II, C II, 
Sill) as well as for deep/self-blended troughs such as O VI, Civ, and Nv. 

For photoionization modeling, actual A^ion] measurements with robust error bars are 
much more useful than lower (or upper) limits. We now demonstrate that the troughs of 
Aim, Pv and Hel* are not heavily saturated and therefore reliable A/" [ion] measurements 
can be extracted for them. In Section [3T21 we show that this is also the case for SlV*. Let 
us first examine the Aim, and Pv doublets, where the oscillator strength (/) of the blue 
transition is twice that of the red transition. From Figure [3] we can infer that in both cases 
the residual intensity of the red doublet component (I r ) is significantly higher than residual 
intensity of the blue doublet component (if,) along essentially the entire absorption trough. 
Such a situation implies the true ionic column density in the trough cannot be much larger 



than 2-3 times that of an AOD estimate (see detailed treatment of this issue in 



Arav et al. 



2005 



20081 ) . Our Table [T] shows that this is indeed the case for Aim, and Pv, where 
we are able to measure robust AT [ion] for these ions using all three absorption models. If 
we define the saturation as S = A^[ion](PL)/A^[ion](AOD) then the level of saturation in 
these troughs (from Tabled]) is less than 1.2 for Aim, and less than 3 for Pv (even when 
taking the large upper error bar on the Pv measurement). Similarly the Hel* is even less 
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saturated (S < 1.1). Under these conditions, the exact behavior of the covering fraction as 
a function of velocity (C(v)) across the trough for the PC model, or that of the exponent 
a(v) of the PL method, do not affect the robustness of the derived A/" [ion]. 

While Si IV does not seem particularly saturated in Figure [3j we have to stress the 
fact that the normalization of the line profile was performed under the assumption that 
the absorbing material covers the emission sources (continuum and emission lines) by 
the same fraction. A closer inspection of the non-normalized O VI and C IV absorption 
troughs, however, reveal the existence of residual emission indicating that at least part 



Arav et al. 


1999a: 


Gabel et al. 


2006) 



by the absorber as observed in other outflows (e.g. 
Failure to account for this observation leads to an underestimation of the total Si IV column 
density since the weaker line emission under the blue absorption trough relative to the red 
would mislead the observer by making the doublet appear non-saturated (see Figure H}. 
For this reason, we report a lower limit on the Si IV column density in Table [TJ Another 
transition that could be affected by this problem is Mgll, although in this case, the weaker 
underlying line emission would only affect the absorption profile of the red component of 
the doublet. Other non-saturated species (Pv, Aim etc.) are located further away from 
significant emission lines in the SDSS J1512+1119 spectrum and are therefore not affected 
by this problem. In the last column of Table [TJ, we report the column densities used in the 
photoionization analysis. Those column densities are selected according to the following 
procedure: we use the values reported in the PC column as the measurements and the PL 
measurement and e rror as the upper err or in order to account for possible inhomogeneities 



in the absorber (see 



Borguet et al. 



20121 ). When only AOD determinations are available, we 



consider the reported values minus the error as a lower limit. 
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Table 1. Computed column densities 

Ion AOD a PC a PL a Adopted b 

10 12 cm" 2 10 12 cm~ 2 10 12 cm- 2 10 12 cm- 2 



Hi 


> 765 




He i* 


687+J 2 


6951}° 


Cn 


214+^ 




Civ 


> 2550 




Nv 


> 3280 




Ovi 


> 3940 




Mgn 


24 8 +0 - 3 


35.6l°j 


Alii 


1 7fi+ - 08 

L - ' °-0.08 




Aim 


42.41°;? 


51.0l°; 8 7 


Sin 


4 05 +a37 




Sim 


no q+0.2 
Z,O.D_q 2 




Si iv 


3911 2 


5291^ 


Pv 


343±| 


4441} 3 


Siv c 







715; 



-27 



51.91};° 



55.71°;? 



> 800 



695 



+211 

-8 



> 765 
695+^ 

> 210 

> 2550 

> 3280 

> 3940 
35.6™ 

> 1.68 



51.01^ 



> 3.7 
> 23.1 
>525 



444 +462 



26600; 



-4000 
-2700 



a) The integrated column densities for the three absorber models. The quote d error arise 
from photon statistics only and are comp uted u sing the technique outlined in iGabel et al. 



(j2005bh . 



b) Adopted values for the photoionization study (see text), 
c) See Section l3~2l for details. 
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3.2. The S IV column density 



The combined good S/N, medium resolution and moderate Ly forest contamination 
in the X-Shooter spectrum of SDSS J1512+1119 allow us to identify and separate the 
absorption t roughs associated w i th Sl V in kinematic component 2 of the outflow. As 



discussed in 



Leighly et al. 



( 2009 



20111 ). S IV is composed of three lines; the ground state 
transition with a wavelength of 1062.66 A, and two transitions arising from an excited state 
(E = 951 cm -1 , referred to as S IV* in the following) at wavelengths 1072.96 and 1073.51 A, 
the excited state being populated at a critical density of 4.7 x 10 4 cm -3 . A useful feature of 
these lines is that the fractional abundance of S IV peaks at a similar ionization parameter 
to that of the ubi quitous Civ, implying that they arise in similar regions of the outflow 



( jDunn et al. 



20121 ) . The prese nce of two excited S IV* transitions with oscillator strengths 



an order of magnitude apart (IHibbert et al. 



20021 ) provides sensitivity to a wide range of 



optical depths before both troughs become fully saturated. 

In Figure HI we present the result of the simultaneous Gaussian fit performed over the 
observed line profiles of the three S IV lines associated with component 2 of the outflow. 
The fit was performed by fixing the wavelength positions of the lines relative to the velocity 
of the component derived from the strong P v A1117.98 line profile as well as imposing an 
identical FWHM for the three Gaussians, leaving only the depth of each component as a 
free parameter. The result of the fit clearly shows a detection of the weak excited A1073.51 
transition, translating to a high S IV* column density. Such a high column density leads 
to the conclusion that the A1072.97 transition is optically thick and that the non-black 
profile observed is due to the partial covering of the emission source by the absorbing 
material. The depth of the resonance A1062.66 line is consistent, within the uncertainties, 
with the depth derived from the A1072.97 excited line implying that the A1062.66 profile 
is also only reflecting a partial covering of the emission source. Note that the presence 
of a trough associated with component 1 of the outflow is observed in SlV around 1055.2 
A and in SlV* A1072.97 around 1065.5 A could affect the result of the fit if the optical 
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Fig. 4. — Result of the simultaneous fit (continuous smooth line) to the three S IV lines 
present in component 2. Within the uncertainties, the fitted profiles of the resonance 
A1062.66 and excited A1072.97 lines are identical, suggesting that the line profiles are satu- 
rated and no column density information can be derived from them. However, the detection 
of the excited 1073.51 A line, a transition ten times weaker than the 1072.97 A line, allows 
us to estimate the real S IV* column density. 
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depth in that system is large enough to produce a significant S IV* A1073.51 trough. 
However, the high covering deduced from other high ionization lines (e.g. C IV, Si iv) 
for that component suggests a low optical depth (r ~ 0.2) associated with the A1072.97 
translating to a r ~ 0.02 for the ten times weaker A1073.51 transition. Such a tiny optical 
depth will not affect the presented modeling of trough 2. Using the Gaussian model of 
the non-saturated A1073.51 line profile along with the fact that the A1072.97 line profile is 
saturated provides an unequivocal determination of the covering fraction across the trough 
(i.e. C = 1 — Ii for saturated lines). We solve the PC model residual intensity equations 
for both S IV* transitions simultaneously and estimate the total S IV* column density to be 
iV[Slv*] ~ 1.6 x 10 16 cm~ 2 . 

In order to determine the total column density in S IV, we also have to be able to 
estimate the column density present in the ground state transition. As detailed above, 
the A1062.66 line is saturated so that no accurate column density can be derived from the 
depth of the line profile. However, knowing the electron density n e of the gas would allow 
us to estimate the total S IV column density by comparing the measured value to models 
predicting the population ratio in the excited state to the ground state as a function of 
n e . The detection of a blend of line s that we identify w ith the Cm* multiplet gives us the 



possibility to do so. As detailed in 



Gabel et al. 



fl2005a[ ). the excited Cm A1175 multiplet 



comprises 6 lines arising from 3 J levels. The J=0 and J=2 levels have significantly lower 
radiative transition probabilities to the ground state than the J=l level an d are thus 
popula ted at much lower densities than the latter. In particular, Figure 5 in 



Gabel et al. 



(I2005al ) shows that the relative populations of the three levels are a sensitive probe to a 
wide range of n e while being insensitive to temperature. In Figure [5j we show several fits of 
the C ill* multiplet assuming a broadened Gaussian optical depth distribution for each line. 
The broadened Gaussian optical depth profile is identical for each line of the blend, only the 
optical depth varyies between the lines. The profile composed of a main central Gaussian 
profile (FWHM = 70 km s _1 ) containing the core of the optical depth (r core ) centered at the 
radial velocity identical to the centroid of the P v All 17.98 feature. Two weaker Gaussians 
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(r = 0.2r core ) with FWHM = 70 km s -1 are added at ±70 km s _1 of the core Gaussian in 
order to produce the broadened wings. The fitting model is tightly constrained by the exact 
kinematic position (derived from the P V trough) and given the kinematic separation of 
the six C ill* components. The fit was repeated for various electron densities in the range 
log(n e ) G [3, 10] with a step of 0.1 dex. The only free parameters for each model were a 
single optical depth and covering fraction. The best fit was found using a x 2 figure of merit 
for a density of log(n e ) ~ 5.4 cm -3 . We estimated conservative error-bars on the electron 
density by noting that for densities lower than log(n e ) ^ 4.8 the population of the J=2 level 
is too low to produce enough absorption (upper left panel of Figure E]) while the absence of 
significant contribution from J=l lines limits the maximum density to log(n e ) ,$,8.1 (lower 
right panel of Figure [5]). Electron number densities within the range 4.8 > log(n e ) > 8.1 
provides acceptable fit to the blend. The density estimated from the C ill* fit is in good 
agreement with the detection of Fe ill* lines from the UV34 multiplet (see Borguet et al. 
2012b, in preparation) which have a critical density of ~ 10 5 cm -3 and with the absence 
of Fe ill* lines from the UV48 multiplet for which the critical density is ~ 10 10 cm -3 
(Bautista 2008, private communication). Assuming the electron density derived above and 
that the level population of S IV* relative to S IV are determined by collisional excitation 
and radiative de-excitation only, we are able to estimate the column density in the S IV 
ground state to be 1.06lo 27 10 16 cm -2 , where the errors reflect the uncertainty on n e . 



4. Photoionization modeling and relative abundance of phosphorus 

4.1. Physical state of the gas assuming solar abundances 

In this section we use the column densities determined for each of the ionic species to 
constrain the physical state (ionization parameter U, and total hydrogen column density 
Njj) of the gas. We model the absorber by a plane parallel slab of gas of constant hydrogen 
number density (uh) and use the spectral synthesis code Cloudy C10. 00 (last described 
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0.2 ; CIII* 1175 multiplet region 
■log(n)=4.5 U 

0.0 



j= 10 12 1 2 




CIII* 1175 multiplet region 1 r 
log(n ) = 5 .4 Best Fit model U 




0.2 ; CIII* 1175 multiplet region 
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0.0 



1 12 1 2 




CIII* 1175 multiplet region r 
log(n)=8.4 U 



1167 1168 1169 1170 1171 1167 1168 1169 1170 1171 
Rest Wavelength (A) Rest Wavelength (A) 



Fig. 5. — Fits to the normalized Cm* multiplet blend (continuous line) for component 2. 
The vertical lines plotted under each sub-component of the multiplet shows the location of 
the centroid of the line and the length of the mark indicates the relative oscillator strength 
of the line. The absorption line located at 1170.6 A is an unrelated line belonging to 
the Ly forest. We present the fit of the blended profile assuming the ratio between the J 
levels for four representative electron number densities. The lower limit on the possible n e 
is constrained by the inability of the model to fit the absorption associated with the longer 
wavelength J=2 transition (upper left). The upper limit on the n e range is fixed by the the 
apparition of significant absorption associated with the J=l level for log(n e ) > 8.1 (lower 
right). The best fit of the blend is obtained for log(n e ) = 5.1 (upper right panel). Other 
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in iFerland et al.l Il998l ) to sol ve the ionizati o n equ ilibrium equations. We assume solar 



abundances for the gas from 



Lodders et al. 



(120091 ) . Due to the lack of constraints on the 



EUV/FUV region of the spectral energy distribution (SED) of SDSS J1512+1119, we adopt 



the "UV-soft" SED proposed for luminous, Radio Quiet quasars in 



Dunn et al. 



Figure 11. The features of t his SED, which is d iscussed in detail in 

I 



20 



Dunn et al 



0): their 



(120 lOh 



departs from the "classical" iMathews fc Ferlandl (119871 ) SED (MF87) by excluding the so 
called UV-bump peaking at FUV energies while keeping an a ox index similar to that of 
MF87. 



We use the grid model approach described in Edmonds et al.l (1201 if ) in order to 
determine the pair of parameters (U,Nh) that best reproduces the observed ionic column 
densities. We only consider column densities determined from non-saturated troughs, 
since the lower limits placed on the column densities from Hi, C IV, N v, O VI are clear 
underestimations of their true values (see Section I3.ip and are consistent with the solution 
derived from the non-saturated lines. We have constraints from 10 ionic species : He I*, 
Cll, Mgn, Al II, Aim, Sin, Sim, Siiv, Pv and SlV. From this set, Pv, SlV, Hel* and 
Aim are the most reliable as their iV[ion] is derived from two or more troughs that are not 
heavily saturated (see Section 13. ip . 

Given that the derived density of the absorbing material is well over the C II critical 
density we conclude that our estimation of the Cll column density is robust within a 
factor of two, since for such a high density, the apparent strength of Cll* A1335.71 is 
approximately twice that of Cll A1334.53, allowing us to derive a PC and PL solution by 
applying a similar treatment to these lines as to resonance doublets. The other lines are 
either singlets for which a robust iV[ion] is difficult to ascertain in principle, or doublets 
that may be more heavily saturated (Siiv, Mgn, see Section \3.1\\ . We finally place an 
upper limit on the column density of the non-detected Fell by scaling the Mgn blue line 
profile template to the 2a noise in the region where the strongest Fell A2382.77 should be 
located and find JV[Feii] < 1.96 x 10 12 cm" 2 . 
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Fig. 6. — Grid model plot of the photoionization modeling of the absorber. Solid lines are 
the adopted iV[ion] measurements from Table [TJ The shaded area represent the measured 
column density of the ion and its errors on these measurements due to the photon noise 
and the systematic uncertainties in the absorber model. Dotted-dashed line represent a 
lower limit on the column density. The black square represents the U and Nh that best 
reproduces the observed column densities (see text). The upper part of the phase plot 
(beyond the hydrogen ionization front characterized by log(iV#) = 23.1 log(LQ) is an artifact 
of our modeling procedure. By coincidence the He I* and S IV lines almost track each other 
between \og(U) = —1.4 and log(£7) = —0.5. We represent the hydrogen ionization front 
(log(JVff) = 23.1 \og(U H )) by a thick black line. 
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Table 2. Photoionization model of the absorber 



Ion log(A obs ) (cm" 2 ) log(A mod ) (cm" 2 ) log(^) 
Adopted a Cloudy 



log (U) ■■■ -0.90 b 



log (N H ) ■ ■ • 21.9 b 



Hi 


> 14.88 


17.51 


2.63 


He I* 




14.86 


0.02 


Cn 


~ 14.33 


14.50 


0.19 


Civ 


> 15.41 


17.91 


2.50 


Nv 


> 15.52 


16.91 


1.40 


Ovi 


> 15.60 


17.81 


2.21 


Mg ii 


io.oo_ 01 


14.34 


0.79 


Alii 


> 12.23 


12.45 


0.22 


Aim 


i o 71 +0.05 
iJ - ' x -0.01 


14.13 


0.42 


Sin 


> 12.57 


12.96 


0.39 


Sim 


> 13.36 


15.09 


1.73 


Si iv 


> 14.71 


16.24 


1.52 


Pv 


14.651^ 


14.65 


0.00 


Siv 


16.431°;™ 


16.45 


0.02 



a) Adopted column densities reported in Table [U Ions with robust measurements are marked in boldface. 

b) Best fit Cloudy model. 
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We present the results for a grid of photoionization models in Figure |6j Visual 
inspection of the figure shows that a minimum ionization parameter of log(£7) > —1.1 is 
required by the Pv constraints. The presence of the low ionization Cll, Sill, and Alii 
requires \og(U) < —0.6: although these low ionization species can also be formed at higher 
U near a hydrogen ionization front in the absorber, such a model would overpredict well 
constrained ions like S IV or P v by a factor > 5. The corresponding total hydrogen column 
density is rather large due to the detection of P v as well as the high column density 
derived from the detailed analysis of the S IV troughs and suggests a thick absorber with 
log(iV#) e [21.5, 21.9] for the range of ionization parameters considered here. Such a high 
column density is also supported by the detection of He I*, which implies an absorber with 
a thickness reaching well into the He II region, while the clear non-detection of the strongest 
Fell lines implies the absence of the formation of a hydrogen ionization front in the slab. 
We determine the best (U, Njj) model by \ 2 minimization of t he difference between the 



Borguet et al 



20121 ) and 



measured column densities and those predicted by Cloudy (see 
find \og(U) = —0.9 and log(iV#) = 21.9. In Table [21 we compare the measured column 
densities to those predicted by the best-fit model. Taking the uncertainties into account, 
the column densities associated with the ions for which we have a robust measurement 
are matched well within a factor of 2 while the Al ill column density is reproduced within 
a factor of 3. The computation shows also an underestimation by factors of hundreds for 
the usual high ionization lines, underlining the fact that AOD measurements yield a poor 
estimate of the actual column densities for these species. 



4.2. Constraining the re 

In this section, we use our knowledge 
measurements of the column densities of P 
abundance of phosphorus in the outflowing 
except phosphorus have solar abundances. 



lative phosphorus abundance 

of the photoionization solution and the 
V and mainly He I* in order to constrain the 
material. Let us first assume that all elements 
Given the photoionization solution determined 
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in Section I4.1[ the abundance of phosphorus relative to helium is constrained by the 
upper-limit on the column density of P v to be ^2 times the solar value. In these models, 
we use a hydrogen number density of 10 4 cm -3 , which is close to the lower limit determined 
in Section 13.21 Increasing the hydrogen number density to 10 8 cm -3 , the upper limit on 
the number density, increases U and Nh by 0.1 dex each and reduces the upper limit on 
the abundance ratio P/P©. We therefore conclude that, for solar abundances of the other 
elements, the phosphorus abundance is approximately solar. 

While the parameter fit using solar abundances is very good, we consider the effects 
of changing metalicity on the phosphorus abundance. For our purposes, we consider a 



gas with metalicity Z/Z© = 1 to have sol ar abundances, noting that the scaling o 



elements with Z is model dependent (e.g. 



Hamann fc Ferland 



1993 



Korista et al. 



heav y 



19961 ). 



Figure [7] is a grid model plot for photoio nization models wit h Z/Z Q « 4. We use the 



abundance scalings provided in Table 2 of 



Ballero et al 



fl2008h for C, N, O, Mg, Si, Ca, and 



Fe. Scalings for the other metals as well as helium are estimated using starburst models in 
Cloudy for Z/Z = 4. For these models, any solution that reasonably reproduces column 
densities of the metals underestimates the column density of helium by a factor ~ 2 for a 
number density of 10 4 cm -3 and less than a factor of 2 for a number density of 10 8 cm~ 3 . 
Increasing the metalicity increases this discrepancy, thus we conclude that Z/Z© ^ 4. We 
find that none of the models accurately reproduce our measured ionic column densities if 
the metal abundances are reduced by a factor ~ 2, implying that the gas is has a metalicity 
Z/Z© ?l 0.5. Constraining the metalicity to 0.5<Z/Z© <4 limits the ionization parameter 
to —1.3 U — 0.5. Comparing phosphorus only to helium, we find that 1/2 <P/P Q < 5. 
The upper limit is overly conservative and is the maximum over- abundance of phosphorus 
that allows the models to produce enough He I*. 

Ionization and thermal structures in the absorber depend on the SED incident on 
the outflowing gas. In the foregoing analysis, we used the UV-soft SED mentioned in 
Section 14.11 We tested several SEDs appropriate for radio-quiet quasars to determine their 
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impact on our main results. The largest changes occurred for SEDs including a substantial 
"UV bump". In particular, using MF87, we find the best fit model with solar abundances 
yields log U = —1.1 and log Nh = 21.6. However, more importantly for our purposes here, 
our diagnostic line ratios change by factors ^ 2 for corresponding models. The major 
change is that the conservative upper limit on P/P© can be as large as 6 for the MF87 
SED. However, we emphasize again that the MF87 SED is not a good approximation to a 
high luminosity radio quiet quasar. 



4.3. The true C IV optical depth 

For outflows that show a significant P v trough, our investigation gives us a unique 
empirical opportunity to contrast the apparent and real optical depth of the C IV trough 
arising from the same outflow. As can be seen in Figure [31 the residual intensity of the 
Civ trough in the deepest part of component 2 is Ij ~ 0.025, which yields an apparent 
optical depth of ln(ij) ~ 4. The real optical depth of Civ at that velocity can be 
estimated using our derived photoionization parameters for the absorber along with the 
knowledge of the column density and abundance of an unsaturated line like P v. Assuming 
solar abundances for simplicity (the abundance of phosphorus does not deviate widely 
from solar, see Section 14.21) we represent in Figure [8] the expected column density ratio 
log(iV[P v]/iV[Clv]) as a function of U and Nh- Using the ionization parameter of 
\og(U) = —0.9 and total hydrogen column density of log(iV#) = 21.9 found in Section fl~Tl 
we find log(A[P v]/N[C iv]) ~ —3.3. The ratio of optical depth for the two transitions is 
then given by: 

TjOvl = iVICl^^v = 10 3,W9 1548^0 ^ 

r[Pv] iV[Pv] (/A) PV 0.45 1117.98 ' v ; 

where / is the oscillator strength and A is the wavelength of the given transition. Using 
the partial covering method on the P v doublet troughs we obtain a real optical depth for 
PvA1117.98 of 2.7. Therefore, equation flU yields a real CivA1548.20 optical depth of 
~3200! Or a value almost 1000 times larger than the apparent optical depth of Civ(~ 4). 
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Fig. 7. — Identical to Figure |6] but assuming a metalicity of Z = 4 (see Section 14.21 for 
details). It is clear that no U and Nh solution can simultaneously satisfy the S IV and He I* 
measurements, as these curves never intersects. 
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Thus the detection of P v troughs will indeed result in heavily saturated C IV line profiles, 
with a real optical depth roughly a thousand times higher than the apparent one. 



5. Summary 



In this paper, we studied the UV outflow associated with the quasar SDSS J1512+1119 
on the basis of new, medium resolution VLT/X-Shooter data. The extended wavelength 
coverage of the instrument allowed us to detect the outflow components in a multitude 
of ionic species. In particular we report the detection of deep Pv absorption troughs in 
kinematic component 2 as well as the detection of Siv and Siv*. A detailed analysis of 
the Siv* line profile allowed us to detect the weak A1073.51 transition revealing a Siv 
column density larger than suggested by the apparent depth of the absorption troughs of 
the A1072.97 transition. 

Photoionization modeling of the absorber revealed that the absorber is thick, though 
the non-detection of significant Fell absorption troughs guarantees the absence of a 
significant amount of H I bound-free opacity. Our accurate determination of the total P v 
S IV He I* and lower ionization species column densities allowed us to characterize the 
physical state of the gas. We find that for the range of ionization parameters relevant for 
the present absorber, the phosphorus abundance relative to helium is consistent with the 
solar value. Using the param eter derived from the photoionization analysis, we show, as 



suggested in 



Hamann 



(119981 ). that a line such as the ubiquitous Civ is heavily saturated. 
The Civ column density derived from the apparent depth of the absorption line profile 
underestimates the column density by a factor of ^1000, providing a very poor estimate of 
its true column density. 

The phosphorus abundance we find is in dis agreement with the extreme phosphorus 



abun dances reported in the early literature (e.g. 



Junkkarinen et al. 



1995 



1997 



Hamann 



19981 ). Other elemental abundances are found to be in agreement with the solar values. 
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-1.2 -1.0 -0.8 -0.6 -0.4 

log(U) 



Fig. 8. — This plot shows contours of the logarithm of the predicted ratio of column density 
between P v and C IV, log(iV[P v]/N[C iv]), as a function of the ionization parameter U and 
the total hydrogen column density of the absorber Njj- Assuming solar abundances and 
the photionization solution found in Section H~Tl (represented by a solid diamond), Cloudy 
models predicts a ratio of column density of log(iV[Pv]/iV[Clv]) ~ —3.3. 



-28 - 



The fact that the abundances are similar to the solar values for an odd (P v) and even 
(S IV") element points to enrichment by relatively "normal" galactic s tellar populatio ns (e.g. 



Hamann 



19971 ) rather than the more exoti c mechanism p roposed by 



would significantly enhance the P/S ratio ( iHamannl 1 1 9981 ) . 



Shieldsl fll996f ) that 
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